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Abstract 

Considering the effective Euler-Heisenberg Lagrangian, i.e., non- linear photon- 
photon interactions, we study the circular polarization of electromagnetic radia- 
tion based on the time-evolution of Stokes parameters. To the leading order, we 
solve the Quantum Boltzmann Equation for the density matrix describing an en- 
semble of photons in the space of energy-momentum and polarization states, and 
calculate the intensity of circular polarizations. Applying these results to a linear 
polarized thermal radiation, we calculate the circular polarization intensity, and 
discuss its possible relevance to the circular polarization intensity of the Cosmic 
Microwave Background radiation. 

PACS: 73.50.Fq, 42.50.Xa, 98.70.Vc 

Introduction. Modern cosmological observations of the cosmic microwave background 
(CMB) radiation provide important evidences to understand our Universe. Cosmolog- 
ical informations encoded in the CMB radiation concerns not only temperature fluc- 
tuations and the spectrum of anisotropy pattern, but also the intensity and spectrum 
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of linear and circular polarizations. It is generally expected that some relevant lin- 
ear and circular polarizations of CMB radiation should be present, and polarization 
fluctuations are smaller than temperature fluctuations [T]. Recently, there are several 
ongoing experiments [2] to attempt to measure CMB polarizations. Theoretical studies 
of CMB polarizations were carried out in Refs. [U H] , and numerical calculations [H E] 
have confirmed that about 10% of the CMB radiation fields are linear polarizations, 
via the Compton and Thompson scatterings of unpolarized photons at the last scatter- 
ing surface (the redshift z ~ 10 3 ). It is important from theoretical points of view to 
understand the generation of CMB linear and circular polarizations. 

In principle, under effects of background fields, particle scatterings and tempera- 
ture fluctuations, linear polarizations of CMB radiation field propagating from the last 
scattering surface can rotate each other and convert to circular polarizations. This is 
described by the formalism of Faraday rotation (FR) and conversion (FC) j7] , and the 
conversion from linear to circular polarizations is given by the time evolution of the 
Stokes parameter V: 

V = 2Uj t (A<f) FC ), (1) 

where A(f)pc is the Faraday conversion phase shift [Sj. 

Refs. [H |9] present the role of background magnetic fields in producing the CMB 
circular polarization, and Atfipc ~ 10~ 19 for micro gauss magnetic fields. In refs. [TU|[TT]. 
the angular power spectrum of CMB circular polarizations and relevant correlations are 
studied in the case that circular polarization is generated by photon-electron scattering 
in the presence of magnetic fields. I is shown that Lorentz symmetry violation in some 
extension of standard model for particle physics 0, d2] and an axion-like cosmological 
pseudoscalar field [T3] can generate the CMB circular polarization. Ref. [H] shows 
that the noncommutative QED with the Seiberg-Witten expansion of fields in the last 
scattering surface can also generate the CMB circular polarization, and the FC phase 
shift A<p FC ~ KT 17 . 

In this letter we show that circular polarizations of radiation fields can be generated 
from the effective Euler-Heisenberg Lagrangian (see review articles in refs. [T ^ [T5 | fT6j). 
By taking into account this effective Lagrangian, and using the Quantum Boltzmann 
Equation [I], we study the time-evolution of the Stokes parameter V. Applying our 
results to the homogeneous CMB radiation with non-vanishing linear polarizations, we 
obtain the upper limit of circular polarization intensity V/Tq < 7 x 10~ 10 in units 
of thermal temperature T of the CMB radiation, and the corresponding A<f>pc < 
1.6 x KT 13 . 
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Stokes parameters. A nearly monochromatic electromagnetic wave propagating in 
the f-direction is described by: 

E x = a x (t) cos[u t - x (t)], E y = a y (t) cos[u t - 6 y (t)], (2) 

where amplitudes a Xsy and phase angles 9 XjV are slowly varying functions with respect to 
the period % = 2ix /ujq. Any correlation between the a x - and a^-components indicates 
polarizations of electromagnetic waves. In a classical description [17], Stokes param- 
eters, which describe polarization states of a nearly monochromatic electromagnetic 
wave, are defined as the following time averages: 

I c = (a 2 x ) + (a 2 y ), 

Qc = (al) - (dy), 

U c = (2a x a y cos(6 x -9 y )), 

V c = (2a x a y sm(9 x - 9 y )), (3) 

where the parameter I c is total intensity, Q c and U c intensities of linear polarizations of 
electromagnetic waves, whereas the V c parameter indicates the difference between left- 
and right- circular polarizations intensities. Linear polarization can also be character- 
ized through a vector of modulus Pl = y/Qc + Uc- ^ i s important to notice that the 
Stocks parameters ([3]) are defined for a monochromatic electromagnetic wave with a 
definite momentum k. Given a linear polarization, one can always transform them to 
a coordinate system where Q c or U c vanishes leaving no circular polarization V c = 0. 
In order to generate a net circular polarization via birefringence there must be some 
special coordinate system so that one linear polarization state propagates differently 
from the other due to interactions. 

In a quantum-mechanical description, Stokes parameters can be equivalently defined 
as follows. An arbitrary polarized state of a photon (|fco| 2 = |k| 2 ), propagating in the 
z-direction, is given by 

|e) = aiexp(z0i)|ei) + a 2 exp(z6» 2 )|e 2 ), (4) 

where linear bases |ei) and |e 2 ) indicate the polarization states in the x- and y-directions. 
Quantum-mechanical operators in this linear bases, corresponding to Stokes parameter, 
are given by 



/ = 


kiXcil 


+ |e 2 > <e 2 1 


Q = 


\ei)(ei\ 


— |e 2 )(e 2 | 


u = 


h)(e 2 | 


+ l e 2)(ei| 



V = z|e 2 )(ei| - i|ei)(e 2 |. (5) 
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An ensemble of photons in a general mixed state is described by a normalized den- 
sity matrix p^ = ( |e») (ej|/trp), and the dimensionless expectation values for Stokes 
parameters are given by 



Q=(Q) 



tr pi = 1, 
trpQ = p n - p 22 , 
trpU = p 12 + P21, 
trpV" = zp 2 i - «P2i, 



(6) 
(7) 
(8) 
(9) 



where "tr" indicates the trace in the space of polarization states. This shows the 
relationship between four Stokes parameters and the 2x2 density matrix p for photon 
polarization states. 

Euler-Heisenberg Lagrangian and circular polarizations. The Euler-Hesinberg 
effective Lagrangian is given as follows: 



£ 



eff 



£n -\- 5 £ , 



(10) 



where the first term £ 
second term S£ 

5£ 



jF^F^ is the classical Maxwell Lagrangian, and the 



a 
90m 4 



(F M ^) 2 + i (F^) 2 



'li: 



where m is the electron mass, F^ u = e^ ua ^F a p (see review articles in refs. [T4"| [T5| IT6]). 

We express the electromagnetic field strength F^ u = d^A v — d u A^ and free gauge 
field in terms of plane wave solutions in the Coulomb gauge |18j . 

d 3 k 



(27r) 3 2fc° 



a t (k)e tll (k)e- lk - x + a\(k)el(k)e 



ik-x 



(12) 



where ej M (/c) are the polarization four-vectors and the index i = 1,2, representing two 
transverse polarizations of a free photon with four-momentum k and k° = |k|. a,i(k) 
[a\(k)} are the creation [annihilation] operators, which satisfy the canonical commuta- 
tion relation 

ai{k),a]{k')] = (27r) 3 2k°8 ij 8 ( - s \k-k'). (13) 

The density operator describing an ensemble of free photons in the space of energy- 
momentum and polarization state is given by 

d 3 p 



P 



tr(p) J (2vr) 3 



Pij(p)a](p)a j (p), 



(14) 
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where Pij(p) is the general density- matrix ([MH]) in the space of polarization states with 



a fixed energy-momentum "p". The number operator D^{k) = a\(k)aj(k) and its 
expectation value is defined by 

(D%{k) ) = tr0£>° (*)] = (27rm0)(2k°)p l3 (k). (15) 

The time evolution of photon polarization states is related to the time evolution of the 
density matrix Pij(k), which is governed by the following Quantum Boltzmann Equation 

(qbe) m, 

(2 7 r) 3 5 3 (0)(2fc )|p u (A ; ) = z([^(t);^.(A ; )])-i J dt( [H°(t); [fl?(0); D% (k)]} >, (16) 

where the interacting Hamiltonian Hj(t) = — 5£ pip . The first term on the right- 
handed side is a forward scattering term, and the second one is a higher order collision 
term. 

It is known that the linear Maxwell Lagrangian £q in Eq. ( flOi) does not gener- 
ate circular polarizations. We attempt to compute the effect of the non-linear Euler 
Heisenberg Lagrangian (TTTT) on the generation of circular polarizations by using QBE 
(I16p . Eq. (TTTj) is perturbatively small, at the order of a 2 , so that we only compute 
the first order of QBE, i.e. the first term in r.h.s. of Eq. (I16j) . and neglect the second 
term which is of the order of a 4 . The contribution from the first term (F^F^) 2 in 
Eq. (TTTT) vanishes, because it is a squared Maxwell action and commutates with the 
number operator . While the second nonlinear term (F^F^) 2 does not commutate 
with the number operator and gives non-vanishing contributions. 

As a result, we approximately obtain the time-evolution equation for the density 
matrix, 

(2 7 r) 3 5 3 (0)2A; ^(A;) « i<[2/J(f), !$(*)]) 



x 



^(2 7 r) 3 5 3 (0)e^e^'A ;7 A ;At ^(A ; )e„,(fc)] 
Pvj {k)8 si - Pis (k)5 l 'i + Psj (k)S Vi - 



AT) 

X / 777-^-^^a^[ e ^(p)^'(p)][^'(p) + Ps'Kp) + SS ( 17 ) 



(27r) 3 2p° J 

The calculations are tedious, but straightforward. We first apply the Wick theorem 
to arrange all creation operators to the left and all annihilation operators to the right, 
then we use the contraction rule 

(4(p>s(p)) = 2p°(2nf5 3 (p-p')p s Ap), (18) 
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to calculate all possible contractions of creation and annihilation operators a\ and aj. 
For example, we calculate the expectation value 

{p\a\ l {p')a s {p)a\ l {q')a l {q) \p) = (jp\a\,{p')a\ l {q')a s {p)ai{q) \p) 

+ 2p°(27r) 3 rt 3 (p - qO<p| 4(p>(<7) \p) 

= 4p°g (27r) 6 5 3 (p - p')5 3 (q - q'W(pW(?) 

+ Vg°(27r) 6 5 3 ( P - q')5 3 (q - p')p s >i(q)[5 s i> + Mp)], 

where the first line results from the Wick theorem and commutation relations (fl3|) . 
while the second line results from all possible contractions ([IB]) of operators a \ and a s . 
Using Eq. (TTTI) . we obtain the time-evolutions for Stocks parameters (EHHJ) as follows: 

i(k) = 0, (19) 
Q(k) = X{[p 21 (k)-p 12 (k)][p 21 (p)-p 12 (p)} 

x [e2u(k)e lf3/ (k)e 1 f } (p)e 2 u'(p) + ei„ (A;) e 2 ^(fc) e 2 p{p)ex uf (p)] 

= -x\v{k)V{p)[e 2y {k)e wl {k)e w {p)e 2yl {p) + e lv (fc)e 2j8 /(A;)e 2j9 (p)e ll/ /(p)]|j(20) 

U(k) = x{[p 21 (k) - p 12 (k))[p n (p) - p 22 {p) ]e 2l/ (k)e 2j3 , (k) e 1/3 (p)e lu t(p)\ 

= x\iV{k)Q(p)e 2u (k)e w ,(k)e w {p)e lv ,(p)^ (21) 

V(k) = X^[p 22 (k) - pu(k)][pi 2 (p) + P2i(p)]eiu(k)€ 1 p l (k)€ 2 p(p)€2u'(p) 

+ [P2i(k) - Pi2(k)][pu(p) - P22(p)}e2u(k)e 2l 3>(k)e ll3 (p)e lul ( y p)} 

= x{Q(k)U(p)e lv {k)e 2pt (k)e 2 pip)e 1 ^(p) 

+ iV{k)Q{p)e2 U {k)e2,y{k)e 1 p{p)e ll// {p)}, (22) 



where k indicates the energy-momentum state of incoming photons in a radiation field, 
and p indicates the energy-momentum states of virtual photons in vacuum, and the 
operator X is defined as following integral overall energy-momentum states p, 

*{-}-^/pfe[^"M MU v] {•■•}■ m 

The I modes represent the ensemble of photons. The Q and U modes represent the 
ensemble of linearly polarized photons, and the V mode represents the ensemble of 
circularly polarized photons. 

Eq. (fTTj) gives an interacting vortex of four photons. Eqs. ( [191122"]) result from the 
tadpole diagram of a photon loop integrating all contributions "p" of virtual photons 
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in vacuum (see Eg. 123]) . This indicates that polarization states of a propagating photon 
with momentum u k" interact with those of virtual photons in vacuum. If the photon 
u k" is not initially polarized, i.e., Q(k) = U{k) = V(k) = 0, then Eqs. (111122]) show that 
the photon u k" propagating through vacuum does not acquire polarizations. Instead, 
the photon u k" is linearly polarized, as if there were a particular orientation of local 
magnetic field. Interacting with this local magnetic field, virtual photons can develop 
circular polarization states, in turn these states back-react with polarization states of 
the photon "/c". As a result, Eqs. ( 1T91122"]) show that the photon "fc" acquires a net 
circular polarization. 

It is important to notice that in the right-handed side of Eq. ( 122]) . the linearly 
polarized modes Q(k) and U(p) are in different momentum states u k" and "p" so that 
they are independent modes. This cannot be made by a coordinate transformation 
so that one of them (U or V) vanishes, as discussed for the Stocks parameters Q or 
(IMS]) for a monochromatic electromagnetic wave with a definite momentum "fc". In 
addition, the right-handed sides of Eqs. (I20|21ll22p show nonlinear interactions between 
Q, U and V modes of a given momentum state "A;" of a radiation field and all possible 
momentum states "p" from vacuum contributions. In these nonlinear interactions, Q, 
U and V modes differently interact with each other leading to circular polarizations, 
i.e., non-vanishing V(k) modes, provided Q(k) and U(p) are non zero. 

In Eq. ( Ti~9]) . 1 = implies in the ensemble of photons, the total intensity of pho- 
tons is constant in time-evolution. In Eqs. (120112 ip and (1221 . the time-evolution Q, 
U and V are given by the combinations of Q, U and V modes, which indicates a 
rotation or conversion between these modes as long as the effective interaction ( ITT]) 
acts. The time-evolution V is proportional to Q and U modes. This indicates that 
an ensemble of linearly polarized photons will acquire circular polarizations due to the 
Euler-Heisenberg Lagrangian (ITT]) . We are interested in considering an initially linear 
polarized electromagnetic radiation propagating through vacuum, and calculating how 
much the intensity of circular polarization can be converted from the intensity of linear 
polarization, due to the non-linear Euler-Heisenberg interaction. 

Intensity of circular polarizations in CMB. We apply our results (I19H22P to the 

ensemble of thermal CMB photons, /bb(p) = l/(e P//T — 1), where the temperature T 
and photon momentum p are in the comoving frame. Thus photon energy and number 
densities are given by: 

f d 3 p 7T 2 4 f d 3 p 2C(3) 3 

£j= J (27F P/bbCp)= 15 ' ™ 7 = J j2^ fBM = (24) 
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and the mean energy for each thermal photon 



&v vr 4 T 



n 7 [30C(3)J 

corresponds to the intensity I c ([3]). We consider that the thermal radiation is initially 
polarized and propagates through vacuum. In order to calculate the final intensity of 
circular polarizations ([22"]) . we approximate 

Q{k)e 2 p'{k)e lv (k) fa C Q f BB (k)S^, U{p)e 2 p{p)e w {p) « C u f BB (p)S (3ul , (26) 

where coefficient Cq(Cu) is the ratio of linear Q(£/)-polarization intensity and total 
intensity. The coefficient Cq represents the linear polarization of the real photon "fc" 
propagating in vacuum Cq < 1, while Cy represents the sum over all contributions of 
linear polarization of virtual photons "p" in vacuum, and Cjj ~ 1. Cq and Cjj are 
independent of each other, because they are associated to different momentum states k 
and p of photons. 

Assuming that the converted intensity of circular polarization is much smaller than 
the intensity of linear polarization, we neglect the second term in the right-handed 
side of Eq. (T2"2"j) . Integrating Eq. (T22J) over all momentum states p [see Eq. (I2"5j) ]. we 
approximately obtain 

V(k) = (^f) (^C v C Q [kf BB {k)i (27) 

where k is the momentum of thermal photons. The form of expression (1271) can be 
understood as the rate of converting the linearly polarized mode to a circularly po- 
larized mode. The linearly polarized Q-mode of momentum state k interacts with 
linearly polarized [/-modes of all momentum states p in vacuum, and converts to the 
circularly polarized V-mode of momentum state k. The factor [kf BB (k)] is due to the 
energy-spectrum of thermal photons and (T/m) 4 comes from the summation over all 
momentum states p ( J2^|) . Integrating Eq. ( 1271) overall energy- momentum states u k" of 
the black-body distribution f BB (k), and normalizing it by the total number- density n 7 
of thermal photons, we obtain 

Finally, multiplying the mean intensity e 7 /n 7 , we obtain 

dV - ^^v^y CuCQ; (29) 



dt V 2 7 Vm7 



n. 



8 



analogously, we redefine the Stokes parameters Q and U in energy units by multiplying 
the mean energy e 7 /n 7 , 




and U =>■ Q, corresponding to those in Eq. (JSJ). 

To estimate the V, we integrate over the comoving time J dt = J dz/H(z), where 
the redshift z e [0,10 3 ], the Hubble function H(z) = H [Q M (z + l) 3 + fi A )] 1/2 for 
fl M ~ 0.3, Qa — 0.7 and H = 75 km/s/Mpc, and the temperature T = T (l + z) 
[T « 2.725fT° = 2.349 x 10~ 4 eV = (0.511cm)" 1 ] in the standard cosmology [191 12Q]. 
From Eqs. (j25|29p we obtain in units of the present CMB temperature T , 

f - -(^)©W° "-* + .m W 

» 7x 10~ w CuC Q < 7x 10" 10 . (31) 

We find that this intensity of circular polarizations is very small, as compared with the 
CMB anisotropy AT /To « 10 -5 [21], which measures the inhomogeneity of the CMB 
radiation. From Eqs. ([I]), ( 1291) and ( 1301) , and we obtain the Faraday conversion phase 
shift 

,22 rp 4 H000 

A0 FC ~ 1.35 (^)p) T C Q / ^(l + 2 ) 5 /i7( 2 ) 



w 1.6 x 10~ 13 C Q < 1.6 x 10" 13 . (32) 

Our results ( )3T|32l) give the upper limit of the effects of non-linear Euler-Heisenberg 
Lagrangian on the intensity of CMB circular polarization. 

Analogously to Ref. [22], which also discusses the effect of photon-photon scatterings 
on the CMB circular polarization, we take the linear polarization coefficient Cq ~ 
3 x 10~ 5 at the last scattering surface, corresponding to ATq/Tq with the maximum 
AT Q k3x lO" 5 ^ in the WMAP data [23], we obtain AV/T w 2.1 x 10" 14 , and 
A(f)pc ~ 4.8 x 10 -18 . This result implies that the Euler-Heisenberg effect on the CMB 
circular polarization could be rather important, compared with those effects mentioned 
in the introductory paragraph. 

Conclusion and remarks. In this letter, by approximately solving the first order of 
Quantum Boltzmann Equation for the density matrix of a photon ensemble, and time- 
evolution of Stokes parameters, we show that propagating photons convert their linear 
polarizations to circular polarizations by the nonlinear Euler-Heisenberg interactions. 
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We discussed this Euler-Heisenberg effect on the circular polarization of CMB photons, 
and showed that this effect is very small, as compared with the present CMB temper- 
ature Tq. Nevertheless, observational studies on such circular polarization are clearly 
warranted. What and how do we need to observe? Is non-zero circular polarization 
at one point in the sky enough, or should there be correlations with the pattern of 
linear polarization? To our knowledge, there is not any current polarization experiment 
that directly measures the CMB circular polarization. However, in the next five years 
considerably more detailed information about the CMB polarization will be delivered 
by the Planck satellite [21] and ground based, high resolution polarization experiments 
such as ACTPol [25], PIXIE [26], SPIDER [27], PolarBear [28], and SPTPol [29]. The 
sensitivity (polarization) ATq/Tq of the Planck satellite for hight and low frequency 
is in the order of 10~ 6 . This seems to be still far from the CMB circular polarization 
AV/Tq rs 2.1 x 10 -14 estimated in this letter. On the other hand, it would be interesting 
to see this Euler-Heisenberg effect on the circular polarization of laser photons [30J . 
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